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ABSTRACT 

We present late-time radio observations of 68 local Type Ibc supernovae, including six events with broad optical 
absorption lines ("hypernovae")- None of these objects exhibit radio emission attributable to off-axis gamma-ray 
burst jets spreading into our line-of-sight. Comparison with our afterglow models reveals the following conclu- 
sions: (1) Less than ~ 10% of Type Ibc supernovae are associated with typical gamma-ray bursts initially directed 
away from our line-of-sight; this places an empirical constraint on the GRB beaming factor of (/)7') ^ 'cor- 
responding to an average jet opening angle, 6 j ^ 0.8 degrees. (2) This holds in particular for the broad-lined 
supernovae (SNe 1997dq, 1997ef, 1998ey, 2002ap, 2002bl and 2003jd) which have been argued to host GRB 
jets. Our observations reveal no evidence for typical (or even sub-energetic) GRBs and rule out the scenario in 
which every broad-lined SN harbors a GRB at the 84% confidence level. Their large photospheric velocities and 
asymmetric ejecta (inferred from spectropolarimetry and nebular spectroscopy) appear to be characteristic of the 
non-relativistic supernova explosion and do not necessarily imply the existence of associated GRB jets. 

Subject headings: gamma-ray bursts — supernovae: specific (SN 2003jd) 



1. INTRODUCTION 

It is now generally accepted that long duration gamma-ray 
bursts (GRBs) give rise to engine-driven relativistic jets as well 
as non-relativistic spherical supernova (SN) explosions. The 
first example of this GRB-SN connection came with the dis- 
covery of the Type Ic supernova, SN 1998bw, associated with 
GRB 980425 {d - 36 Mpc; Galama et al. 1998; Pian et al. 
2000). The unusually fast photospheric velocities and excep- 
tionally bright radio emission of SN 1998bw indicated ~ 10^^ 
erg of kinetic energy and mildly relativistic ejecta (bulk Lorentz 
factor, r - 3; Kulkarni et al. 1998; Iwamoto et al. 1998; 
Li & Chevaher 1999; Woosley, Eastman & Schmidt 1999). 
In comparison with other core-collapse events {Eke ~ 10^' 
erg and ejecta speeds, v <; 0.1c), SN 1998bw was consid- 
ered a hyper-energetic supernova ("hypernova"; Iwamoto et 
al. 1998). Broad optical absorption lines were also observed 
in the Type Ic SNe 2003dh and 20031w, associated with the 
cosmological GRBs 030329 and 031203, indicative of com- 
parably large photospheric velocities (Matheson et al. 2003b; 
Malesani et al. 2004). Together, these observations appear to 
suggest that broad spectral features are characteristic of GRB- 
associated SNe. 

In addition to events with prompt gamma-ray emission, the 
GRB-SN connection also implies the existence of "orphan" su- 
pernovae whose relativistic jets are initially beamed away from 
our line of sight (Rhoads 1999; Paczynski 2001). Since the 
discovery of SN 1998bw, several broad-lined SNe have been 
identified locally {d < 100 Mpc) and are currently estimated to 
represent ~ 5% of the Type Ibc supernova (SNe Ibc) population 



(Podsiadlowski et al. 2004). Given their spectral similarity to 
the GRB-associated SNe, it has been argued that local broad- 
lined supernovae can be used as signposts for GRBs. 

Thus, associations with poorly-localized BATSE bursts have 
been invoked for the broad-lined SNe 1997cy, 1997ef and 
1999E'' (Germany et al. 2000; Turatto et al. 2000; Wang & 
Wheeler 1998; MazzaH, Iwamoto & Nomoto 2000; Rigon et 
al. 2003). In addition, association with off-axis GRBs have 
also been claimed. In the case of SN 2002ap, broad optical ab- 
sorption lines and evidence for mildly asymmetric ejecta (based 
on spectropolarimetry measurements) were interpreted to sup- 
port an off-axis GRB jet (Kawabata et al. 2002; Totani 2003, 
but see Leonard et al. 2002). 

More recently, an off-axis GRB model has been proposed 
for SN2003jd, for which photospheric velocities upward of 
40,000 km s"' were measured at early time (Filippenko, Foley 
& Swift 2003; Matheson et al. 2003a). More intriguingly, late- 
time {t ^ 400 days) spectra showed double-peaked emission 
lines of light-elements, attributed to an asymmetric explosion 
(Kawabata et al. 2004). Mazzali et al. (2005) ai-gue that these 
observations can be understood if SN 2003jd was accompanied 
by a highly coUimated GRB jet initially directed ^ 70 degrees 
away from our line-of-sight. 

Regardless of viewing angle, however, strong afterglow 
emission eventually becomes visible as the decelerating GRB 
jets spread laterally and the emission becomes effectively 
isotropic. As the jets spread into our line-of-sight, a rapid 
increase of broadband synchrotron emission is observed on a 
timescale of a few weeks to several years. This late-time emis- 
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sion is most easily detected at long wavelengths (Perna & Loeb 
1998; Levinson et al. 2002; Waxman 2004). Targeting local 
Type Ibc supernovae with late-time radio observations has thus 
become the preferred method to search for evidence of off-axis 
GRBs (Stockdale et al. 2003; Soderberg, Frail & Wieringa 
2004). 

Using early radio observations {t <^ 100 days) we have al- 
ready Umited the fraction of SNe Ibc harboring on-axis (or 
mildly off-axis) GRBs to be ^ 3% (Berger et al. 2003a). In this 
paper, we present late-time (f ^ 0.5 to 20 yr) radio observations 
for 68 local Type Ibc supernovae, including SN 2003jd and five 
additional broad-lined events, making this the most comprehen- 
sive study of late-time radio emission from SNe Ibc. We use 
these data to constrain the SN fraction associated with ORB 
jets regardless of viewing angle assumptions, constraining even 
those initially beamed perpendicular to our line-of-sight. 

2. RADIO OBSERVATIONS 

2.1. TypeIcSN2003jd 

SN2003jd was discovered on 2003 October 25.2 UT within 
host galaxy MCG -01-59-021 (di ^ 81 Mpc; Burket et al. 
2003). In Table 1 we summarize our radio observations for 
SN2003jd, spanning 8-569 days after the explosion^. All ob- 
servations were conducted with the Very Large Array^ (VLA) 
in the standard continuum mode with a bandwidth of 2 x 50 
MHz centered at 4.86, 8.46 or 22.5 GHz. We used 3C48 and 
3C147 (J0137H-331 and J0542H-498) for flux calibration, while 
J2323-032 was used to monitor the phase. Data were reduced 
using standard packages within the Astronomical Image Pro- 
cessing System (AIPS). 

No radio emission was detected at the optical SN position 
during our early observations. Our radio limits imply that 
SN 2003jd was a factor of ^ 100 less luminous than SN 1998bw 
on a comparable timescale. We conclude that SN2003jd, like 
the majority of SNe Ibc, did not produce relativistic ejecta along 
our line-of-sight. 

We re-observed SN 2003jd at ? ~ 1 .6 yrs to search for radio 
emission from an off-axis GRB jet. No emission was detected, 
implying a limit of F^, < 45 /uJy {3a) at 8.46 GHz. 

2.2. Late-time data on Local Type Ibc Supernovae 

We supplement these data with late-time {t ^ 0.5-20 year) 
radio observations for 67 local (di <^ 200 Mpc) SNe Ibc, sum- 
marized in Table 2. Eleven objects were observed at moder- 
ately late-time as part of our on-going VLA program to char- 
acterize the early {t <^ 100 days) radio emission from SNe Ibc 
(Soderberg et al, in prep). The remaining 54 objects were ob- 
served on a later timescale (f ^ 1 year) and were taken from 
the VLA archive^. We note that five of these supernovae (SNe 
1997dq, 1997ef, 1998ey, 2002ap, 2002bl) were spectroscopi- 
cally observed to have broad optical absorption lines, similar to 
SN 1998bw. 

All VLA observations were conducted at 8.46 GHz (ex- 
cept for SN1991D at 4.86 GHz) in the standard continuum 
mode with a bandwidth of 2 x 50 MHz. Data were reduced 
using AIPS, and the resulting flux density measurements for 
this sample of SNe Ibc is given in Table 2. With the excep- 

'Here we assume an approximate explosion date of 2003 October 21 UT, based 
^The Very Large Array is a facility of the National Science Foundation operated 
'http://e2e.nrao.edu/archive/ 



tion of SN 2001 em, from which radio emission from the non- 
relativistic, spherical supernova ejecta is still detected at late- 
time (Stockdale et al. 2005; Bietenholz & Bartel 2005, but see 
Granot & Ramirez-Ruiz 2004), none of the SNe Ibc show ra- 
dio emission above our average detection limit of ^ 0.15 mJy 
(3(7). In comparison with SN 1998bw, only SN2001em shows 
a comparable radio luminosity on this timescale. These results 
are consistent with the earUer report by Stockdale et al. (2003). 

In Figure 1 we plot the radio observations for this sample 
of SNe Ibc, in addition to late-time radio data for SN1954A 
(Eck, Cowan & Branch 2002) and SN 1984L (Soderberg, Frail 
& Wieringa 2004). 



3. OFF-AXIS MODELS FOR GAMMA-RAY BURSTS 

3.1. An Analytic Approach 

Waxman (2004) present an analytic model for the late-time 
radio emission from a typical GRB viewed significantly away 
from the collimation axis. In this model, the GRB jet is ini- 
tially characterized by a narrow opening angle, Oj ~ few de- 
grees, while the viewing angle is assumed to be large, 6'obs ^ 1 
radian. As the jet sweeps up circumstellar material (GSM) and 
decelerates, it eventually undergoes a dynamical transition to 
sub-relativistic expansion (Frail, Waxman & Kulkarni 2000). 
The timescale for this non-relativistic transition is estimated at 
tNR ~ 0.2(£'5i/no)'^^ yr (~ 0.3£'5i/A* yr) in the case of a homo- 
geneous (wind- stratified) medium (Waxman 2004). Here, £51 
is the beaming-corrected ejecta energy normalized to 10^' erg 
and «() is Ihc circumstellar density of the homogeneous medium 
(interstellar medium; ISM) normaUzed to 1 particle cm"-*. For 
a wind-stratified medium, A* defines the circumstellar density 
in terms of the progenitor mass loss rate, M, and wind velocity, 
v„, such thatM/47rvH, = 5 x 10"A* g cm"', and thus A* = 1 for 
M = 10"^ Mq yr ' and v,„ = 10^^ km s"' (Li & ChevaUer 1999). 

Once sub-relativistic, the jets spread sideways, rapidly in- 
tersecting our line-of-sight as the ejecta approach spherical 
symmetry. At this point the afterglow emission is effec- 
tively isotropic and appears similar to both on-axis and off- 
axis observers. The broadband emission observed from the 
sub-relativistic ejecta is described by a standard synchrotron 
spectrum, characterized by three break frequencies: the syn- 
chrotron self-absorption frequency, the characteristic syn- 
chrotron frequency, Vm, and the synchrotron cooling frequency, 
Uc. On timescales comparable to the non-relativistic transition, 
i^a and iy„ are typically below the radio band while is gener- 
ally near the optical (Frail, Waxman & Kulkarni 2000; Berger, 
Kulkarni & Frail 2004; Frail et al. 2005). Making the usual 
assumption that the kinetic energy is partitioned between rela- 
tivistic electrons and magnetic fields (e^ and cb, respectively), 
and that these fractions are constant throughout the evolution 
of the jet, Waxman (2004) estimate the radio luminosity of the 
sub-relativistic, isotropic emission to be 

on pre-explosion images (Burket et al. 2003). 
under cooperative agreement by Associated Universities, Inc. 
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i„ „ 8.0xlO«(i^)(ii)'„fE„ (1) 

VlOGHzy V^NR/ 
for the ISM case, while for a wind- stratified medium 



VlOGHz/ V^NR/ 

Here it is assumed that the electrons are accelerated into a 
power-law distribution, A^(7) oc 7"'' with p = 2.0. These equa- 
tions reveal that the strength of the non-relativistic emission is 
strongly dependent on the density of the circumsteUar medium 
(especially in the case of a wind) and is best probed at low fre- 
quencies. 

While this analytic model provides robust predictions for the 
afterglow emission att > ?nr. it does not describe the early evo- 
lution or the transition from relativistic to sub-relativistic ex- 
pansion. At early time, the observed emission from an off-axis 
GRB is strongly dependent on the viewing angle and dynamics 
of the jet. To investigate this early afterglow evolution and the 
transition to sub-relativistic expansion, we developed a detailed 
semi-analytic model, described below. 

3.2. A Semi-analytic Model 

In modehng the afterglow emission from an off-axis GRB 
jet, we adopt the standard framework for a adiabatic blastwave 
expanding into either a uniform or wind stratified medium (Sari 
1997; Granot & Sari 2002). We assume a uniform, sharp-edged 
jet such that Lorentz factor and energy are constant over the 
jet surface. The hydrodynamic evolution of the jet is fully de- 
scribed in Oren, Nakar & Piran (2004). As the bulk Lorentz 
factor of the ejecta approaches F ^ 1, the jets begin to spread 
laterally at the sound speed Our off-axis light-curves are ob- 
tained by integrating the afterglow emission over equal arrival 
time surface. We note that these resulting hght-curves are in 
broad agreement with Model 2 of Granot et al. (2002) and are 
consistent with Waxman's analytic model (§3.1) on timescales, 

t ^ fNR- 

Over-plotted in Figure 1 are our off-axis models calculated 
for both wind-stratified and homogeneous media at an observ- 
ing frequency of i^obs = 8.46 GHz. We assume standard GRB 
parameters of £51 = A* = n = 1, = = 0.1, p = 2.2 and 
9j = 5°, consistent with the typical values inferred from broad- 
band modeling of GRBs (Panaitescu & Kumar 2002; Yost et 
al. 2003; Chevalier, Li & Fransson 2004). We compute model 
light-curves for off-axis viewing angles between 30 and 90 
degrees. As clearly shown in the figure, the majority of our 
late-time SNe Ibc hmits are significantly fainter than all of the 
model Ught-curves, constraining even the extreme case where 
^obs=90°. 

'"Since the spreading behavior of relativistic GRB jets is poorly constrained 
conservative assumption since it produces the faintest off-axis hght-curves. 



4. SN2003JD: CONSTRAINTS ON THE OFF-AXIS lET 

Based on the double-peaked profiles observed for the nebular 
lines of neutral oxygen and magnesium, Mazzali et al. (2005) 
argue that SN 2003jd was an aspherical, axisymmetric explo- 
sion viewed near the equatorial plane. They suggest that this 
asynmietry may be explained if the SN explosion was accom- 
panied by a tightiy coUimated and relativistic GRB jet, initially 
directed ^ 70 degrees from our line-of-sight. This hypothesis is 
consistent with the observed lack of prompt gamma-ray emis- 
sion (Hurley et al. 2003) as well as the absence of strong radio 
and X-ray emission at early time (Soderberg, Kulkami & Frail 
2003; Watson ef a/. 2003). 

Our radio observation of SN 2003jd at f ^ 1 .6 years imposes 
strong constraints on the putative off-axis GRB jet. While the 
early data constrain only mildly off-axis jets (^obs ^ 30°), our 
late-time epoch constrains even those jets initially directed per- 
pendicular to our line-of-sight. As shown in Figure 1, our radio 
limit is a factor of ^ 200 (^ 20) fainter than that predicted for 
a typical GRB expanding into a homogeneous (wind- stratified) 
medium, even in the extreme case where ^obs ~ 90°. Given the 
assumption of typical GRB parameters, we conclude that our 
late-time radio limit is inconsistent with the presence of an off- 
axis GRB jet. We note that the model assumptions and physical 
parameters of our off-axis afterglow light-curves are identical 
to those adopted by Mazzali et al. (2005). 

We next explore the range of parameters ruled out by our 
deep radio limits. As shown in Equations 1 and 2, the luminos- 
ity of the late-time emission is a function of the ejecta energy, 
the density of the circumsteUar medium and the equipartition 
fractions. To investigate the effect of energy and density on 
the late-time radio luminosity, we fix the equipartition fractions 
to €e = eB = 0.1, chosen to be consistent with the values typi- 
cally inferred from afterglow modeling of cosmological GRBs 
(Panaitescu & Kumar 2002; Yost et al. 2003). 

In Figure 2, we illustrate how each radio epoch for SN 2003jd 
maps to a curve within the two-dimensional parameter space of 
kinetic energy and circumsteUar density for an off-axis GRB. 
Here we adopt our semi-analytic model (§3.2) for a wind- 
stratified medium, along with a typical electron index of p = 2.2 
and a viewing angle of 6'obs = 90°; the faintest model for a 
given set of equipartition fractions. By comparing the luminos- 
ity Umit for SN 2003jd at a particular epoch with the off-axis 
model prediction for that time, we exclude the region of pa- 
rameter space rightward of the curve since this region produces 
a jet which is brighter than the observed limit. The union of 
these regions represents the total parameter space ruled out for 
an associated GRB. As shown in this figure, the total excluded 
parameter space extends from A* ^ 0.03 and E ~ lO''^ to 10^^. 

We compare these constraints with the beaming-corrected ki- 
netic energies and CSM densities for 18 cosmological GRBs 
(Table 3). Here we make the rough approximation that A* w no; 
a reasonable assumption for circumsteUar radii near lO'^ cm. 
As shown in Figure 2, these GRBs span the region of parameter 
space roughly bracketed byA, -0.002 to 100 and £ - 2 x 10'*'' 
to 4 X 10^^ The majority of the bursts (13 out of 18) fall within 
the excluded region of parameter space for SN 2003jd. We con- 
clude that SN2003jd was not likely associated with a typical 
GRB at a confidence level of ^ 72%. 

by observations, we assume negligible spreading during this phase. We adopt this 
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5. LOCAL TYPE IBC SUPERNOVAE: FURTHER CONSTRAINTS 

While physical parameters atypical of the cosmological 
GRB population can be invoked to hide an off-axis GRB for 
SN 2003jd, it is exceedingly unlikely for atypical parameters to 
dominate a large statistical sample of SNe Ibc. Motivated thus, 
we searched for off-axis GRBs in the 67 local Type Ibc SNe 
for which we have compiled late- time ~ 0.5 - 30 yr) radio 
observations. Applying the method described in §4 we produce 
exclusion regions in the £51 -A* parameter space for each SN. 
Figure 3 shows the resulting contours for all 68 SNe, includ- 
ing SN 2003jd and five broad-hned events. For the twenty SNe 
with early radio limits (Berger et al. 2003a; Berger, Kulkami 
& Chevalier 2002) we combine late- and early-time data to pro- 
vide further constraints. 

In Figure 4 we compile all 68 exclusion regions to quantify 
the £51 - A* parameter space constrained by this statistical sam- 
ple. Contours map the regions excluded by incremental frac- 
tions of our sample. As in the case of SN2003jd, all curves 
rule out bursts with A, > 1 and E > 10^" erg. Moreover, 50% 

exclude A* 0.1 and E ;> lO"*** erg. For comparison, the mean 
ejecta energy and CSM density values for cosmological GRBs 
are £ « 4.4 x 10^° erg and A* = no « 1-2. 

Focusing on the subsample of broad-hned SNe, we empha- 
size that our deep limits rule out both putative GRB jets directed 
along our line-of-sight (e.g. SN 1997ef) as well as those which 
are initially beamed off-axis (e.g. SN2002ap and SN2003jd). 
In particular, the large exclusion region for SN 2002ap (see Fig- 
ure 3) implies that an extremely low CSM density, less than 
A* ^ 3 X 10"^, is needed to suppress the emission from an asso- 
ciated GRB. This is a factor of ^ 10 below the density inferred 
from modeling of the early radio emission (Berger, Kulkarni 
& Chevaher 2002) and we therefore conclude that an off-axis 
GRB model is inconsistent with our late-time observations of 
SN2002ap. In Figure 4 we show that this entire sample of six 
broad-hned SNe rule out bursts with energies £ ^ 10"*^ erg, and 
50% even rule out E ~ lO'*^ erg (all assuming a typical A* = 1). 

We next address the limits on an association with GRBs de- 
fined by the cosmological sample (Table 3). For each SN in our 
sample we calculate the fraction of observed GRBs that lie in its 
exclusion region. We then determine the probability of finding 
null-detections for our entire sample by calculating the product 
of the individual probabilities. We find that the probability that 
every Type Ibc supernova has an associated GRB is 1.1 x 10"'°. 
We further rule out a scenario in which one in ten SNe Ibc is 
associated with a GRB at a confidence level of ~ 90%. For 
the broad-lined events alone we rule out the scenario that every 
event is associated with a GRB at a confidence level of ^ 84%. 
Confidence levels are shown as a function of GRB/SN fraction 
in Figure 5. 

6. DISCUSSION AND CONCLUSIONS 

We present late-time radio observations for 68 local Type 
Ibc supemovae, including six broad-lined SNe ("hypemovae"), 
making this the most comprehensive study of late-time radio 
emission from SNe Ibc. None of these objects show evidence 
for bright, late-time radio emission that could be attributed to 
off-axis jets coming into our line-of-sight. Comparison with 
our most conservative off-axis GRB afterglow models reveals 
the following conclusions: 

(1) Less than 10% of Type Ibc supemovae are associated 



with GRBs. These data impose an empirical constraint on the 
GRB beaming factor, (/^' ), where fh = (1 -cos 9j). Assuming a 
local GRB rate of ~ 0.5 Gpc'^ yr"' (Schmidt 2001; Perna, Sari 
& Frail 2003; Guetta, Piran & Waxman 2005) and an observed 
SNe Ibc rate of - 4.8 x 10"* Gpc'^ yr"' (Marzke et al. 1998; 
Cappellaro, Evans & Turatto 1999; Folkes et al. 1999), we con- 
strain the GRB beaming factor to be (/^') ^ x lO"'. Adopting 

a lower limit of (/^'} > 13 (Levinson et al. 2002), the beam- 
ing factor is now observationally bound by (/^') « [13- W], 
consistent with the observed distribution of jet opening angles 
(Frail et al. 2001; Guetta, Piran & Waxman 2005). 

(2) Despite predictions that most or all broad-lined SNe Ibc 
harbor GRB jets (Podsiadlowski et al. 2004), our radio ob- 
servations for six broad- lined events (SNe 1997dq, 1997ef, 
1998ey, 2002ap, 2002bl and 2003jd) reveal no evidence for as- 
sociation with typical (or even sub-energetic) GRBs. While un- 
usual physical parameters can suppress the radio emission from 
off-axis jets in any one SN, it is unlikely that all six broad-lined 
events host atypical GRBs. We observationally rule out the sce- 
nario in which every broad-lined SN harbors GRB jets with a 
confidence level of ^ 84%. 

(3) While low CSM densities (e.g. A* <^ 0.1) can suppress 
the emission from off-axis GRB jets, such values are inconsis- 
tent with the mass loss rates measured from local Wolf-Rayet 
stars (0.6-9.5 x 10"^ M© yr"'; Cappa, Goss & van der Hucht 
2004), thought to be the progenitors of long-duration gamma- 
ray bursts. 

(4) While we have so far considered only the signature from 
a highly coUimated GRB jet, these late- time radio data also 
impose constraints on the presence of broader jets and/or jet 
cocoons. As demonslralcd by GRBs 980425 and 030329, the 
fraction of energy coupled to mildly relativistic and mildly col- 
Umated ejecta can dominate the total relativistic energy budget 
(Kulkarni et al. 1998; Berger et al. 2003b). Less sensitive to 
to the effects of beaming and viewing geometry, broad jets are 
more easily probed at early time (f 100 days) when the emis- 
sion is brightest. Still, we note that the majority of our late- 
time radio limits are significantly fainter than GRBs 980425 
and 030329 on a comparable timescale, thus constraining even 
mildly relativistic ejecta. 

These conclusions, taken together with the broad spectral 
features observed for GRB-associated SNe 1998bw, 2003dh 
and 20031w, motivate the question: what is the connection be- 
tween GRBs and local Type Ibc supemovae? While current 
optical data suggest that all GRB-SNe are broad-lined, our late- 
time radio observations clearly show that the inverse is not 
trae: broad optical absorption fines do not serve as a rehable 
proxy for relativistic ejecta. This suggests that their observed 
large photospheric velocities and asymmetric ejecta are often 
merely characteristics of the non-relativistic SN explosion and 
thus manifestations of the diversity within SNe Ibc. 
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Table 1 

Radio Observations of SN2003jd 



Date Obs 




GHz 


Fv,SA6 GHz 


Fv,12.5 GHz 


(UT) 


(days) 


(MJy) 


(MJy) 


(MJy) 


2003 Oct 29 


8 


±52 


±34 


±58 


2003 Nov 4 


14 




±77 




2003 Nov 15 


25 




±74 




2005 May 12 


569 




±15 





^Assuming an explosion date of 2003 October 21 UT, based on 
pre-explosion images (Burket et al. 2003). 

''AH flux densities are given as Ict (rms). 
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Table 2 

Late-time Radio Observations of Type Ibc Supernovae 



SN name 


Host Galaxy 


Distance" 


Explosion Date 


lAUC 


Date Obs 


At 


Flux density 






(Mpc) 


(UT) 


(#) 


(UT) 




(days) 


(/uJy) 


1983N 


NGC 5236 


7.2 


1983 Jun 26 


3835 


2003 Oct 


17 


7416 


±124 


1985F 


NGC 4618 


1.1 


1985 May 14 


4042 


2003 Oct 


17 


6730 


±37 


1987M 


NGC 2715 


18.9 


1987 Aug 31 


4451 


2003 Oct 


17 


5891 


±34 


1990B 


NGC 4568 


32.0 


1989 Dec 23-1990 Jan 20 


4949 


2003 Oct 


17 


5032 


±24 


1990U 


NGC 7479 


33.7 


1990 Jun 28- Jul 27 


5063 


2003 Oct 


17 


4845 


±44 


1991A 


IC 2973 


45.5 


1990 Dec 6-13 


5178 


2003 Oct 


17 


4695 


±35 


1991D 


Anon. 


173 


1991 Jan 16 


5153 


1992 Oct 


11 


633 


±49 


1991N 


NGC 3310 


14.0 


1991 Feb 20-Mar29 


5227 


2003 Oct 


17 


4604 


±59 


1991ar 


IC 49 


65.0 


1991 Aug 3 


5334 


2003 Oct 


17 


4458 


±30 


19941 


NGC 5194 


6.5 


1994 Apr 1-2 


5961 


2003 Oct 


17 


3486 


±42 


1994ai 


NGC 908 


21.3 


1994 Dec 8-16 


6120 


2003 Oct 


17 


3336 


±40 


1996D 


NGC 1614 


68.1 


1996 Jan 28 


6317 


2003 Oct 


15 


2818 


±153 


1996N 


NGC 1398 


19.7 


1996 Feb 16-Mar 13 


6351 


2003 Oct 


15 


2787 


±36 


1996aq 


NGC 5584 


23.2 


1996 Jul 30 


6465 


2003 Oct 


17 


2635 


±32 


1997B 


IC 438 


44.3 


1996 Dec 14 


6535 


2003 Oct 


15 


2498 


±33 


1997C 


NGC 3160 


99.2 


1996 Dec 18-1997 Jan 14 


6536 


2003 Oct 


17 


2481 


±37 


1997X 


NGC 4691 


15.7 


1997 Jan 16-Feb 2 


6552 


2003 Oct 


17 


2457 


±22 


1997dc 


NGC 7678 


49.6 


1997 Jul 22 


6715 


2003 Oct 


17 


2278 


±36 


1997dq* 


NGC 3810 


14.0 


1997 Oct 16 


6770 


2003 Oct 


17 


2192 


±30 


1997ef* 


UGC 4107 


49.8 


1997 Nov 16-26 


6778 


2003 Oct 


17 


2156 


±33 


1998T 


NGC 3690 


44.3 


1998 Feb 8-Mar 3 


6830 


2003 Oct 


17 


2066 


±482 


1998ey* 


NGC 7080 


69.0 


1998 Nov 1-Dec5 


6830 


2003 Oct 


17 


1794 


±32 


1999bc 


UGC 4433 


90.2 


1999 Jan 30 


7133 


2003 Oct 


15 


1721 


±25 


1999di 


NGC 776 


70.2 


1999 Jul 6 


7234 


2003 Oct 


17 


1564 


±35 


1999dn 


NGC 7714 


39.7 


1999 Aug 10-20 


7241 


2003 Oct 


1 / 


1524 


±42 


1999ec 


NGC 2207 


38.9 


1999 Aug 24 


7268 


ZUUj (JCl 


1 c 

i J 


1515 


±45 


1999eh 


NGC 2770 


27.6 


1999 Jul 26 


7282 




1 7 


1544 


±19 


2000C1' 


NGC 2415 


53.8 


1999 Dec 30-2000 Jan 4 


7348 




1 7 


1385 


±35 


2000F 


IC 302 


84.5 


1999 Dec 30-2000 Jan 10 


7353 


9003 Oct 




1383 


±26 


2000S 


MCG -01-27-2 


85.8 


1999 Oct 9 


7384 


9003 Oct 


17 


1469 


±30 


2000de 


NGC 4384 


35.6 


2000 Jul 13 


7478 


2003 Oct 


17 


1191 


±40*= 


2000ds 


NGC 2768 


19.4 


2000 May 28 


7507 


2003 Oct 


17 


1237 


±34 


2000dv 


UGC 4671 


57.7 


2000 Jul 4 


7510 


2003 Oct 


17 


1200 


±32 


2001Bt 


IC 391 


22.0 


2000 Dec 25-2001 Jan 4 


7555 


2003 Oct 


15 


1021 


±30 


200 IMt 


NGC 3240 


50.6 


2001 Jan 3-21 


7568 


2003 Oct 


17 


1008 


±54 


2001bbt 


IC 4319 


66.3 


2001 Apr 15-29 


7614 


2003 Oct 


17 


908 


±35 


200 Ich 


MCG -01-54-1 


41.6 


2001 Mar 23 


7637 


2003 Oct 


17 


938 


±32 


2001ci'f' 


NGC 3079 


15.8 


2001 Apr 17-25 


7638 


2003 Oct 


17 


909 


±188 


2001efl^ 


IC 381 


35.1 


2001 Aug 29-Sep 9 


7710 


2003 Oct 


17 


774 


±33 


2001ejt 


UGC 3829 


57.4 


2001 Aug 30 


7719 


2003 Oct 


15 


778 


±36 


zUUlem 




OJ.O 


/UUl oep lU-lJ 


/ /ZZ 


2003 Oct 


17 




yU/ ± JO 


200 list 


NGC 1961 


56.0 


2001 Dec 14-23 


7782 


2003 Oct 


15 


668 


±31 


2002jt 


NGC 3464 


53.0 


2002 Jan 15-21 


7800 


2003 Oct 


17 


637 


±23 


2002ap^t 


NGC 628 


9.3 


2002 Jan 28-29 


7810 


2003 Oct 


17 


626 


±29 


2002bl*t 


UGC 5499 


67.8 


2002 Jan 31 


7845 


2003 Oct 


17 


624 


±39 


2002bmt 


MCG -01-32-1 


78.0 


2002 Jan 16-Mar 6 


7845 


2003 Oct 


17 


630 


±39 


2002cp1' 


NGC 3074 


73.4 


2002 Apr 11-28 


7887 


2003 Oct 


17 


547 


±29 


2002hft 


MCG -05-3-20 


80.2 


2002 Oct 22-29 


8004 


2003 Oct 


17 


357 


±48 


2002hot 


NGC 4210 


38.8 


2002 Sep 24 


8011 


2003 Oct 


17 


389 


±35 


2002hyt 


NGC 3464 


53.0 


2002 Oct 13-Nov 12 


8016 


2003 Oct 


17 


354 


±23 



Table 2 — Continued 



SNname 


Host Galaxy 


Distance'' 


Explosion Date'' 


lAUC 


Date Obs 




Flux density'' 






(Mpc) 


(UT) 


(#) 


(UT) 


(days) 


(/«Jy) 


2002hzt 


UGC 12044 


77.8 


2002 Nov 2-12 


8017 


2003 Oct 17 


344 


±29 


2002jit 


NGC 3655 


20.8 


2002 Oct 20 


8025 


2003 Oct 17 


362 


±18 


2002jjt 


IC340 


60.1 


2002 Oct 22 


8026 


2003 Oct 15 


360 


±30 


2002jpt 


NGC 3313 


52.7 


2002 Oct 20 


8031 


2003 Oct 17 


362 


±37 


2003A 


UGC 5904 


94.4 


2002 Nov 23-2003 Jan 5 


8041 


2003 Jun 15 


182 


±69 


20031 


IC2481 


76.0 


2002 Dec 14 


8046 


2003 Jun 15 


183 


±60 


2003aa 


NGC 3367 


43.1 


2003 Jan 24-31 


8063 


2003 Jun 15 


138 


±70 


2003bm 


UGC 4226 


113.7 


2003 Feb 20-Mar 3 


8086 


2003 Jun 15 


109 


±57 


2003bu 


NGC 5393 


86.1 


2003 Feb 10-Marll 


8092 


2003 Jun 15 


110 


±71 


2003cr 


UGC 9639 


155.2 


2003Feb24-Mar31 


8103 


2003 Nov 2 


233 


±59 


2003dg 


UGC 6934 


79.1 


2003 Mar 24- Apr 8 


8113 


2003 Nov 2 


215 


±90 


2003dr 


NGC 5714 


31.7 


2003 Apr 8-12 


8117 


2003 Nov 2 


206 


±58 


2003ds 


NGC 3191 


132.9 


2003 Mar25-Apr 14 


8120 


2003 Nov 2 


212 


±67 


2003el 


NGC 5000 


80.2 


2003 Apr 19-May 22 


8135 


2003 Nov 2 


180 


±65 


2003ev 


Anonymous 


103.3 


2003 Apr 9-Jun 1 


8140 


2003 Nov 2 


180 


±61 


2003jd*t 


MCG -01-59-21 


80.8 


2003 Oct 16-25 


8232 


2005 May 12 


569 


±15 



^Distances have been calculated from the host galaxy redshift. We assume Ho = 72 km s ' Mpc ' , Q,m = 0.27 and Q,a = 0.73 . 

''Explosion dates are given as a range constrained by the date of discovery and the most recent pre-discovery image in 
which the SN is not detected. In the cases where pre-discovery images are not available (or not constraining), we adopt the 
spectroscopically derived age estimate with respect to maximum light and assume a typical light-curve rise time of ~ 21 days 
to provide a rough estimate of the explosion date. 

'^Calculated using the average of the explosion date range in Column 4. 

''ah observations were conducted at 8.46 GHz (except for SN 1991D at 4.86 GHz) and uncertainties represent \a rms 
noise. 

^We detect a ~ 0.7 mJy source located ~ 2 arcsec from the optical position for SN 2000de. Due to the significant positional 
offset, we assume the radio source is not associated with the SN. 

*Broad absorption Unes observed spectroscopically. 
I^Early radio limits reported in Berger et al. (2003a) 

*Early radio data for SN 2002ap was reported by Berger, Kulkami & Chevalier (2002) 
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Table 3 

Physical Parameters for GRBs 



GRB 


z 


Ejecta Energy^ 
(xlO^i erg) 


Density'' 
(A* = n) 


Density Profile 


Reference'^ 


970508 


0.835 


? 4+1-4 


2 4+1-4 


ISM 


1 


980329 


^3.9 


1 1+0.26 
^•^-0.46 


20+1" 
^^-10 


ISM 


2 


980425 


0.0085 


0.012 


0.04 


Wind 


3 


980519 


It 






ISM 


4 


980703 


0.966 


o c+1.26 
J.J_0.42 




ISM 


2 


990123 


1.60 


1 5+°-^3 


o.ooi9!!!:«!»^ 


ISM 


4 


990510 


1.619 


140-49 




ISM 


4 


991208 


0.706 




18+22 


ISM 


4 


991216 


1.02 


11+°' 




ISM 


4 


000301c 


2.03 


31+0 03 
"••'■'-o.os 


27!^ 


ISM 


4 


000418 


1.118 


3.4 


0.02 


ISM 


5 


000926 


2.066 


9 n+0 34 
^■"-0.2 


161 


ISM 


2 


010222 


1.477 


0.51 


1.7 


ISM 


4 


011121 


0.36 


0.2 


0.015 


Wind 


6 


020405 


0.69 


0.3 


^0.07 


Wind 


6 


020903 


0.251 


0.4 


100 


ISM 


7 


030329 


0.168 


0.67 


3.0 


ISM 


8 


031203 


0.105 


0.017 


0.6 


ISM 


9 



'^Corrected for beaming. Errors bracket the 90% confidence interval. 
''Errors bracket the 90% confidence interval. 

'^References: 1. Berger, Kulkarni & Frail (2004); 2. Yost et al. (2003); 3. Li & 
Chevalier (1999); 4. Panaitescu & Kumar (2002); 5. Berger et al. (2001); 6. Cheva- 
her, Li & Fransson (2004); 7. Soderberg et al. (2004a); 8. Berger et al. (2003b); 9. 

Soderberg ef fl/. (2004b) 



i Redshift Unknown; z=l was assumed. 



11 



10 



32 



N 



m 
O 



^1 



T3 

cd 



1024 



"1 1 1 — I — I I I I I 1 1 1 — I — I I I I I 1 1 1 — I — I I I I 

Wind 

ISM 



T Type Ibc SN 
• Broad-lined SN 
(§) SN2003jd 




_L/ I I I I M 



J I I I I I I 



J I I I I I I I 



10 102 103 104 

Time (days) 

Fig. 1.— Late-time radio limits (Scr) are shown for 62 local Type Ibc supernovae (grey triangles) and six broad-lined SNe (cir- 
cles/arrows) including SN2003jd (encircled dot). For SN2001em, we adopt the spherical SN emission as an upper limit on the 
emission from an off-axis GRB jet. Afterglow models for a typical GRB (£51 = A* = no = 1, = €b = 0.1, Oj = 5°, p = 2.2) are shown 
for both wind-stratified (solid lines) and homogeneous (dashed lines) media at viewing angles of 6'obs = 30° (thickest), 60° (thick) and 
90° (thin) away from the initial coUimation axis of the jet. As discussed in §3.2, we make the conservative assumption that lateral jet 
spreading does not begin until after the non-relativistic transition, at which point the jets spread sideways at the sound speed. All of 
the broad-lined SNe, including SN 2003jd, are significantly fainter than the model light-curves, regardless of the CSM density profile 
or viewing angle. 
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Density (A* = n) 

Fig. 2. — Constraints on an off-axis GRB associated with SN 2003jd. Each of the four 8.46 GHz observations is mapped to a contour 
(thin black Hnes) in the two-dimensional parameter space of kinetic energy and CSM density for an off-axis GRB. As discussed 
in §4, we adopt a conservative off-axis model (wind- stratified CSM, Ce = cb = 0.1, p = 2.2, ^obs = 90°) which produces the faintest 
light-curves in Figure 1. We exclude the grey hatched region of parameter space rightward of each contour since this region produces 
a jet which is brighter than the observed 3(7 limit. The union of these regions (thick black line) represents the total parameter space 
ruled out for an associated GRB. For comparison, we also show are the energy and density values for 17 GRBs (filled circles) from 
Table 3, inferred from broadband modeUng of the afterglow emission. For the GRBs without reported error estimates, we adopt 
a factor of two for the uncertainty of both parameters. We also include the parameters for GRB 980425 (unfilled circle) derived 
from the early radio emission emitted from the quasi-spherical SN component, but note that the lack of strong radio emission from 
SN 2003jd at early time indicates that it is not similar to GRB 980425. 
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Density (A* = n) 



Fig. 3. — Late-time radio limits for 62 SNe Ibc (grey contours) and six broad-lined SNe (black contours, labels) are used to constrain 
the parameters (£51 and A*) of putative off-axis GRBs as discussed in §5. For SN 2001em, we adopt the spherical SN emission as an 
upper limit on the emission from an off-axis GRB jet. As discussed in §4, we adopt a conservative off-axis model (wind- stratified 
CSM, ec = eB = 0.l, p = 2.2, 6'obs = 90°) which produces the faintest light-curves in Figure 1. The region rightward of each SN contour 
is ruled out since it implies an off-axis GRB which is brighter than the observed 3 a hmit. For the 21 SNe with additional early- time 
radio observations, we determine the union of parameter space which is excluded (see Figure 2). Energy and density values for 17 
cosmological GRBs (filled circles) and GRB 980425 (unfilled circle) are shown for comparison. For the GRBs without reported 
error estimates, we adopt a factor of two for the uncertainty of both parameters. It is clear that most of the GRBs lie in the excluded 
regions for each SN, implying that if these events are associated with off-axis GRBs, their densities are unusually low. 
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Fig. 4. — Top panel: Regions of £51 - A* parameter space for off-axis GRBs ruled out by the SNe Ibc sample in Figure 3. Con- 
tours/shading depict the regions ruled out by fractions of the sample, labeled in increments of 10%. The darkest shading corresponds 
to the region excluded by 100% of the sample. The statistical mean £51 and A* values for cosmological GRBs is shown as a barred 
point. Bottom panel: Same as above, except only including the broad-lined SNe. 
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Fig. 5.— Using late-time radio data for 68 local SNe Ibc (Table 2) and the physical parameters inferred for cosmological GRBs 
(Table 3) we constrain the fraction of SNe harboring GRBs. We find that we can rule out the the scenario in which every supernova 
(circles) has an associated ORB (GRBs/SNe=l) at a confidence level of ~ 100% and can rule out a fraction of ^ 0. 1 with a confidence 
level of ~ 90%. For the broad-lined supemovae (stars) we rule out the scenario in which every event harbors a GRB with a confidence 
level of ~ 84%. 



